Adiabatic invariants for the dynamics of rotating 3-dimensional fluid 
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The 3-dimensional fluid dynamics under rapid rotation with beta effect is shown to possess three 
adiabatic-type invariants; their presence leading to the energy accumulation in zonal jets. 
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A. Introduction 

The present note is an addendum to the recent pa- 
per which described three adiabatic-type invariants 
for the rotating shallow water dynamics; their presence 
implying the energy transfer towards zonal jets. The 
addendum is due to the remarks of L. Smith }2]. She 
remarked that it would be interesting if the extra con- 
servation held for 3-dimensional (3D) fluid dynamics and 
expressed belief that the 3D extension would be possi- 
ble. However, this seemed highly improbable, because 
the adiabatic conservation — related to the conserva- 
tion in triad resonances — is essentially tied to the two- 
dimensionality. And even if there were 3D extension of 
some adiabatic invariant, say enstrophy, the resulting in- 
variant would hardly be positive definite; and so, it would 
be useless for the inverse cascade argument. 

Separately, application of [1] to the real oceans and 
atmospheres of rotating planets encounters the follow- 
ing natural question. In a fluid layer of finite depth H, 
the shallow water motions [which are 2-dimensional (2D)] 
can generate — due to nonlincarity — the 3D inertia 
waves. One can argue that the length of a typical Rossby 
wave significantly exceeds the depth H, and so, the en- 
ergy exchange between Rossby waves and inertia waves 
can be disregarded due to the disparity of scales. But 
what about short Rossby waves? Their wavelength can 
be comparable to the depth H. Therefore, it appears 
crucial whether the extra conservations can be extended 
to the 3D fluid dynamics. 

Surprisingly, the 3D extension has turned out to be 
true. It is due to the fact discovered by H.P. Greenspan 
long time ago [3j. In the absence of the beta effect (i.e. 
on the /-plane), he showed that a 3D fluid flow can be 
decomposed into a single 2D geostrophic mode and an 
infinite number of 3D inertia waves [see Eq. ^ below] 
and found that "the inertia modes do not produce triad- 
resonance response on the geostrophic mode" . It is im- 
portant that there is only one geostrophic mode; oth- 
erwise, the enstrophy (which is, roughly speaking, the 
quadratic zonal momentum) would not be positive defi- 
nite, and the inverse cascade argument would fail. The 
third adiabatic conservation would also fail, and so would 
do the balance argument for the energy transfer towards 
zonal jets. 

Consider the dynamics of 3D inviscid fluid of constant 



density; the fluid velocity v(ir, y, z, t) satisfying the equa- 
tions 

v t + (v • V)v + f x v = -VII, (la) 
V • v = 0. (lb) 

Euler's equation (pj,) involves Coriolis force; f is the dou- 
ble angular velocity; II is the fluid pressure (centrifu- 
gal force included) divided by the density. The fluid is 
assumed to occupy a layer between parallel horizontal 
planes, rotating about vertical axis z. The boundary 
condition is vanishing of the vertical fluid velocity at the 
planes. 

The system ([1]) is known to conserve the energy and he- 
licity. Here we show that in the presence of the beta effect 
— when f = (/o + (3oy)i slowly depends on the longitu- 
dinal coordinate y (fg is a reference value of the Coriolis 
parameter) — there are three integrals that are conserved 
by the Euler dynamics ([TJ adiabatically. [The adiabatic 
conservation means approximate conservation over long 
time. It is often due to the slow time- dependence of some 
parameters, but here — like in [l| — it is due to the small 
nonlinearity, compared to the beta-effect, and the slow 
spatial dependence f(y).] 

The three adiabatic invariants are described in the next 
Section [Bj and then in Section [C] it is noted that their 
presence implies the emergence of zonal jets. 



B. Adiabatic invariants 

First, assume that the fluid equations are written in di- 
mensionless form with time scale 1/ fo, length scale L and 
velocity scale foL. Then the dimensionless typical wave 
amplitude is the Rossby number Uq/LJq, with Uq being 
typical dimensional velocity. Dimensionless Coriolis pa- 
rameter is f = (1 + /3y)z with dimensionless parameter 
(3 = L/V/o- 

The general solution of the linearized (about v = 0) 
system (JXJ) with (3 — is a superposition of normal modes, 
of which there are two types |3J: (i) an infinite number 
of inertial oscillations (with amplitudes a m and real non- 
zero frequencies Lo m ) and (ii) a single geostrophic mode 
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U = (U, V), which is two-dimensional and steady 

oo 

v(x, y, z, t) = U(ar, y) + a m V m (x, y, z)e"™' (2) 

771 — — OO 

[the sum is over m = ±1,±2, . . . (m = is excluded); 
0J-m = —Mm, v m = V* m , a m = a* m ]. When the non- 
linearity is present (but small) and (3 ^ (but small), 
the amplitudes a m and the geostrophic mode U acquire 
slow time-dependence. And the Euler system (JT|) can 
be written as a system of evolution equations for U and 
amplitudes a m - 

For the nonlinear system with (5 — 0, "the inertia 
modes do not produce triad-resonance response on the 
geostrophic mode" U; in other words, the coupling co- 
efficient between the geostrophic mode (which has zero 
frequency) and two inertia modes (with frequencies to m 
and uj n ) vanishes whenever u> m + uj n = 0. 

When /3 7^ 0, the geostrophic mode becomes the quasi- 
geostrophic mode, or Rossby wave with dispersion law 

fi(k) = -^ (k 2 =p 2 +q 2 ) 

[k = (p, q) is the horizontal wave vector] . As well, the 
inertia wave frequencies acquire some small corrections 
(corresponding to the beta correction in the Coriolis pa- 
rameter). Now solving for a m (t) and U(a;, y, t) by pertur- 
bations, we encounter small denominators, which are not 
exactly canceled by the corresponding numerators (and 
could produce secular terms). 

Let us estimate the effect of nonzero (3 on the small 
denominators. We will consider the mid-latitudinal beta- 
plane, when /3o//o ~ V-^Oj Rq being the Earth radius. 
There are two types of triad interactions between Rossby 
and inertia waves: (i) two Rossby waves and one inertia 
wave, and (ii) two inertia waves and one Rossby wave. 

The first type interactions are clearly far from be- 
ing resonant, and so, do not produce small denomina- 
tors. Indeed, if H <C Rg, then the Rossby frequencies 

~ fioL are much smaller than the inertia frequen- 
cies oj m — f k z / y^k 2 + fej, with vertical wave number 
k z ~ mn/H (m = 1, 2, . . .) and horizontal wave number 
k ~ tt/L. 

The second type can be resonant and produce small de- 
nominators. However, the parameter (3 ^ will change 
small denominators only slightly: the Rossby frequen- 
cies and the change in the inertia frequencies are both 
~ PoL <C /o (at mid-latitudes). Therefore, it is possible 
to add 0(/3)-corrections to the coupling coefficients to 
have exact cancellation of the small denominators, and 
then show that the compensation for these corrections 
remains within the accuracy of adiabatic conservation. 
So, the rapidly rotating 3D fluid dynamics adiabatically 
conserves the invariants of the 2D shallow water dynam- 
ics. 

The present work is somewhat parallel to the study 
[l[ (containing many details) of adiabatic conservation 
in the rotating shallow water dynamics. It also has two 



types of waves: Rossby waves and inertia-gravity waves; 
the latter also do not produce triad-resonance response 
on the former. For the shallow water case, it was possi- 
ble to consider situation when the correction {3qL to the 
Coriolis parameter exceeds /o (which happens, e.g., near 
the equator). In this situation, parameter (3 should be 
small compared to y 'gH /L 2 (g is the gravity constant). 
It is unclear at present whether such situation admits 
extra conservation for the rotating 3D fluid dynamics. 
The adiabatic invariants are integrals 

I=\J X k \( k \ 2 dk, (3) 

where £k is the 2D Fourier transform of the 2D vorticity 
£(x, y,t) =V X — U y , which can be found by averaging the 
vertical component of the 3D vorticity V x v along the 
vertical direction 

I f H 

C(x,y,t) = — J (v x ~u y )dz. 

The kernel is an arbitrary linear combination of the 
following 3 functions 

X k = i (4a) 

X k = 1, (4b) 

p 2 ( P 2 + 5g 2 ) 
Xk = jfcio ' ( 4c ) 

The function (j4^,) determines the energy of the quasi- 
geostrophic mode, its enstrophy corresponds to the func- 
tion (0b). The function (HJ;) appeared in [1] for the dy- 
namics described by a single quasigeostrophic equation 
in the limit when the Rossby radius of deformation p ap- 
proaches infinity. For the shallow water system, the limit 
p — > oo leads to small denominators of two kinds; both 
can be canceled [l|; and the same function (Ufc) deter- 
mines the extra invariant in the rotating shallow water 
dynamics. 

The adiabatic conservation of the integral ([3]) means 
that variation AI = I(t) — 1(0) remains small (compared 
to I) for a long time t of many wave periods. Of course, 
the adiabatic conservation holds if / is sufficiently small. 
However, A J can be small even if / is large: I(t) can 
oscillate in time, as usually the case for adiabatic con- 
servation in the theory of dynamical systems; and this is 
also the case here. 



C. Emergence of zonal jets 

In rapidly rotating 3D fluids, the energy is transferred 
towards 2D flow. This was demonstrated in numeri- 
cal simulations @; this paper describes several possible 
mechanisms of such transfer beyond triad resonances. So, 
the energy accumulates in the 2D quasigeostrophic mode. 

Any of the kernels (HJ) is positive, and the correspond- 
ing invariant Q is positive definite. The invariants Q 
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with kernels (j4^,b) imply the inverse cascade of the energy 
of the quasigeostrophic mode. The presence of all three 
adiabatic invariants (jl^ibc) implies [l| that the energy of 
the quasigeostrophic mode flows to large scale zonal flow 
(not just to large scales, but to the region of the k-plane 
with polar angles 9 = arctan(g/p) close to ±7r/2). 

Thus, the rapidly rotating 3D dynamics transfers en- 
ergy towards zonal jets. 

D. Conclusion 

We have found three adiabatic invariants of the rapidly 
rotating 3D fluid dynamics (namely the energy, enstro- 
phy, and extra invariant of the quasigeostrophic mode) . 
Their presence together with the two-dimensionalization 
(see Sec. [C]) implies the transfer of energy towards zonal 
jets. 

It would be interesting to consider adiabatic invari- 



ants for 3D fluid layer with a free surface. This situation 
presents a combination of the present work with [l[ and 
would lead to the general Rossby wave dispersion, with 
a finite Rossby radius of deformation p. In that general 
situation, the 3D fluid dynamics still seems to have the 
three adiabatic invariants, which would still lead to the 
energy transfer towards zonal jets — except for the long 
wave limit, when the forcing scale (at which the energy 
is supplied) is much bigger than p (the relation between 
L and p can be arbitrary). 
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